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Abstract 

Building on the Renormalization Group (RG) method the beam-beam interaction 
in circular colliders is studied. A regularized symplectic RG beam-beam map, 
that describes successfully the long-time asymptotic behavior of the original sys- 
tem has been obtained. The integral of motion possessed by the regularized RG 
map has been used to construct the invariant phase space density (stationary 
distribution function), and a coupled set of nonlinear integral equations for the 
distributions of the two colliding beams has been derived. 
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1 Introduction 



The problem of coherent beam-beam interaction in storage ring colliders is one of the most 
important, and at the same time one of the most difficult problems in contemporary acceler- 
ator physics. Its importance lies in the fact that beam-beam interaction is the basic factor, 
limiting the luminosity of a circular collider. Nevertheless, some progress in the analytical 
treatment of the coherent beam-beam interaction has been made - Qj, it is still far from 
being completely understood. In most of the references available the basic trend of analy- 
sis follows the perturbative solution of the Vlasov-Poisson equations, where the linearized 
system is cast in the form of an eigenvalue problem for the eigenmodes. 

An important question, which still remains unanswered is how to determine the invariant 
phase space density (equilibrium distribution function) if such exist. One possible way to 
approach this problem is to find an integral of motion (at least approximately) under certain 
conditions. Then the invariant density can be expressed as a generic function of the integral 
of motion. An attempt in this direction has been made by Alexahin || , who used the Deprit 
algorithm to determine the integral of motion (new action variable). 

In the present paper we develop a novel approach to the beam-beam interaction in circular 
colliders, based on the Renormalization Group (RG) method ||. Originally this method 
has been proposed as a singular perturbation technique for differential equations. Naive 
perturbation expansions |7|] are well-known to produce secular terms, thus limiting the range 
of validity of the perturbation solution. The basic idea of the RG method is to remove secular 
or divergent terms by renormalizing the integration constants of the lowest order perturbative 
solution. Its extension to discrete symplectic maps is however not straightforward, and 
should be performed with care. Here we follow the regularization procedure outlined in 
the paper by Goto and Nozaki ||. As shown in || the naive RG map, obtained as a 
result of renormalization of the lowest order solution preserves the symplectic symmetry only 
approximately, and does not describe the long-term behavior of the original map correctly. 
The symplecticity is recovered by a process of "exponentiation", yielding a symplectic RG 
map together with an explicit expression for the nonlinear tune shift. An alternative version 
of the RG method, based on the envelope technique |§ has been applied to study non 
symplectic maps. 

The paper is organized as follows. In the next Section we derive the one-dimensional 
nonlinear beam-beam map. In Section 3 the regularized RG map and its integral of motion 
are obtained. The integral of motion thus found is further used in Section 4 to derive a set 
of coupled integral equations of Haissinski type for the invariant phase space density. 

2 The Nonlinear Beam-Beam Map 

We begin with the one-dimensional model of coherent beam-beam interaction in the vertical 
(q) direction, described by the Hamiltonian 
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where the normalized beam-beam potential V k (q; 9) satisfies the Poisson equation 



d 2 V k 
dq 2 

and 



Airjdpf3- k {q,p;e), (2.2) 



Afc _ Rr e N 3 _ k 0* kq 1 + Ao/3( 3 - fc )o _ 2Rr e N 3 _ k (3* kq 

Here, (k — 1,2) labels the counter- propagating beams, is the azimuthal angle, Xk = Rfikq 
is the derivative of the phase advance with respect to 6, R is the mean machine radius, r e 
is the classical electron radius, iV lj2 is the total number of particles in either beam, f31 q is 
the vertical beta-function at the interaction point, and L kx is the horizontal dimension of 
the beam ribbon. In addition, the distribution function fk(q,P] 6) is a solution to the Vlasov 
equation 

dfk | ^ Kp d fk dH k df k _ Q 
86 dq dq dp 

In order to build the iterative beam-beam map we formally solve the Hamilton's equations 
of motion 

dq dp 

where the prime implies differentiation with respect to the spatial variable q. By defining 
the state vector 



q 
p 



(2.6) 



we can rewrite Eq. ( |2.5| ) in a vector form 

z = K(#)z + F(z;#), (2.7) 

where 

R M=(4o')' (2 ' 8) 

Performing a linear transformation defined as 

z = M{9)1 (2.9) 
where the matrix is a solution of the linear equation with a supplementary initial condition 

M = K(9)M, M(9 )=I, (2.10) 
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we write the equation for the transformed state vector £ as follows: 

'^=M- 1 (9)F(z;9), f(0 o )=z„. 
Equation ( |2.11|) can be solved directly to give 



:2.iD 



z(0) = M(9)z + dTM(9)M- L (r)F(z(ry,r). 



(2.12) 



It can be easily checked that the matrix of fundamental solutions M. to the unperturbed 
problem is of the form 



M{9) 



cos [Xk{9) - Xfc(0o)] s in \Xk{9) - Xk%)\ 
- sin [xk(9) - Xk(0o)] cos [xk(9) - Xk(9 )) 



(2-13) 



so that 



z(0) = M(9)z - X k drS p (T)Vl(q(r); r 



00 



sin \xk{9) -X&(r)] 
cos [xk(0) - Xk(r)] 



(2.14) 



Applying the above expression (|2.14 ) in a small e-interval 9 & (9q — e, 9 + e) around the 
interaction point (located at 9q) and then taking the limit e^Owe obtain the kick map 



Qkick — Qo, 



Pkick = Po~ AfcV^(g ) 



(2.15) 



In order to obtain the rotation map in between successive kicks we apply once again expres- 
sion qn3D in the interval 9 <E (9 + e, 9 + 2n - e) 



g rot = ?kick cos 2nv k + p kick sin 2^, 

Prot = -gkick Sin 27TI/ fc + Pkick COS 27TI/ fc . 



(2.16) 



Combining Eqs. ( |2.15| ) and ( |2. 16| ) we finally arrive at the one-turn beam-beam map 



Qn+i = Qn cos2vrz/ fe + \p n - X k V k \q n )} sm2nv k , 
Pn+i = -q n sm2TTiy k + [p n - X k V k (q n )] cos 2^- 



(2.17) 



It is important to note that the one-turn beam-beam map (|2.17|) is symplectic, since its 
Jacobian determinant is equal to unity 



det 



d{q n+ i,p n+1 ] 
d(q n ,p n ) 



det 



cos 2nis k — X k V k sin 27rz^ sin 2i^v k 
— sin 2nis k — X k V k cos 2Tru k cos 2nv k 



1. 



:2.i8i 
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3 Renormalization Group Reduction of the Beam- 
Beam Map 

The one-turn beam-beam map, derived in the previous Section can be further simplified by 
eliminating the canonical momentum variable p from ( |2.17| ). Multiplying the first of Eqs. 
( 2.1 7| ) by cos 2-KP k , multiplying the second one by — sin 27rz^., and summing the two equations 
up we find 



q n+ i cosu k - p n+1 smu k = q n , 

where 

u k = 2irv k . 

Using Eq. fl3.1p we obtain a second order difference equation 

Cq n = q n+ i - 2q n cosuJ k + g n -i = -t\ k V k \q n ) smu k , 



(3.i; 



(3.2) 



(3.3) 



where e is a formal small parameter (set to unity at the end of the calculations), taking into 
account the fact that the beam-beam kick is small and can be treated as perturbation. 



Next we consider an asymptotic solution of the map ( |3.3j ) for small e by means of the 
RG method. The naive perturbation expansion 



(3.4) 

when substituted into Eq. ( |3.3|) yields the perturbation equations order by order 

£q<® = 0, (3.5) 



Cq^ = -X k V k ^)smuj k , 
Cqi? = -\ k £ ) V k "(q^)smuj k , 



CS = -A* 



V k "'(q^)+q^V k "(q^ 



sin uj k , 



Solving Eq. ( |3.5| ) for the zeroth order contribution we obtain the obvious result 

A k e iuJkTl + c.c. = 2\A k \ cos (u k n + (f> k ), 



JO) 
Sin 



f = iA k e iUkn + c.c. = -2\A k \ sin (u k n + <j> k ), 



(3.6) 
(3.7) 

(3.8) 

(3.9) 
(3.10) 



where A k is a complex integration constant, whose amplitude and phase are \A k \ and cp k 
respectively. 
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Let us assume for the time being that the beam-beam potential Vk(q) is a known function 
of the vertical displacement q. In what follows it will prove efficient to take into account 
the fact that the beam-beam potential Vk(q) is an even function of the coordinate q. Odd 
multipole contributions to Vk(q) will give rise to a shift in the closed orbit, and can be easily 
incorporated in the the calculations presented below. It is straightforward to check that the 
Fourier image of the beam-beam potential Vk(X), defined as 

oo oo 

V k (q) = — J dXV k (X)e iX «, V k (X) = J dqV k {q)e- iX \ (3.11) 

— oo — oo 

retains the symmetry properties of V k (q), that is: 

V k (-X) = V k (X). (3.12) 



Using the expansion |fL0| , [TT| 



e 



iz cos ip 



E i m J m (z)e mip , (3.13) 



where J m {z) is the Bessel function of the first kind of order m, and the explicit form of the 
zero order solution ( pT9| ) we find 

oo 

V k \q^) = E C ( k M) Al M -V^ M -^ + c.c. (3.14) 

M=l 

Here the coefficients are functions of the amplitude \A k \ and are given by the expression 
C { k M \\A k \) = - \ > /dAA^(A) l 7 2M _ 1 (2A|^|). (3.15) 

Similarly for the second derivative of the beam-beam potential V k " \Qnj> entering the second 
order perturbation equation ( |3.7| ) we have 



oo 

VZ(q<?) = + E V^Al M e^ + c.c, (3.16) 



M=l 



where 



oo 

V ( k \\A k \) = — fdXX 2 V k (X)J (2X\A k \), (3.17) 

7T J 



vl M \\A k \) = - { -^~ fdXX 2 V k (X)J 2M (2X\A k \). (3.18) 

7T \A k \ J 



A/+1 00 
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From the recursion property of Bessel functions jTOL III 



Ju-l(z) + J v+1 {z) = —J u (z) 

z 



(3.19) 



we deduce an important relation to be used later 



V W_ V (N + l) lAkl 2 = {2N+1)c (N + l)^ 



(3.20) 



The solutions of the perturbation equations ( |3.6| ) and (|3.7|), taking into account ( |3.20| ) are 
given by 



In 



(1) _ l^k 



—C^A k e iuJkn + Sm Uk V C^ M+1) ■ A 2 k M+l e i{ - 2M+l)LJkn + 
2 2 M=l 



c.c. 



(3.21) 



\ 2 ■ \ 2 oo 

8 v y 4 iV=l 



+ 



N=l 



n + % 



sin(2iV + l)u k 



cosu k — cos(2iV + l)u> k 



A2N+1J(2N+I)u> k n 



AfcSin uJk (Q) 



(AT+l) 



j coscufc — cos(2iV + l)u k 



2 „• 2 , oo oo 



+ 



Aj,sin 



E E 



^M+l)p(A0 



I j^ijv=i COSUJ k - cos[2(M + N) + l]u k 



A 2(M+N)+l e i[ 2 (M+N)+l] 



Lu k n 



+ 



XI sin 2 oo k 00 ' 00 ' 



E E 



^(M+lWAO | , ,4iV 
<-fc ^fe l^fcl 



4 m^in=i coso; fc -cos[2(M- N) + lju^ 



A 



2{M-N)+1 i[2(M-N)+l]u> h n 



+ c.c. (3.22) 



where the summation in the last term of Eq. ( |3.22| ) is performed for M ^ N, and 



C 



(N+l) 



c 



(N+l) 



cosu k — cos(2A r + l)to k 



(3.23) 



To remove secular terms, proportional to n and n 2 we define the renormalization trans- 
formation A k — > A k (n) by collecting all terms proportional to the fundamental harmonic 

e iw k n 



A k (n) = A k + e^fc^A k 
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~-C k ^ 2 (n 2 + in cot u k 

H x 



i\\ sin u k x 



AT=1 



(3.24) 



Solving perturbatively Eq. (|3.24|) for in terms of A k {ri) we obtain 



A fc fn) 



A discrete version of the RG equation can be defined by considering the difference 

iX 



A k (n + l)-A k (n) = e^C { k 1) A k 



(3.25) 



\2 

k -Cl 1)2 (2n + l + icotu k ) 



i\\ sin ui k 



JV=1 



(JV+lWiV), A ,47V 
U k n k 



A k . 



(3.26) 



Substituting the expression for A k in terms of A k (n) [see Eq. ( |3.25|) 1 into the above Eq. 
(|3.26| ) we can eliminate the secular terms up to 0(e 2 ). The result is 



A k {n + l) 



1 + e— Q. ; — e — C k ' (1 + ^cotwfcj 



2 A| sin cu fc 



iV=l 





4JV 


A fc (n) 





A fc (n) 



(3.27) 



This naive RG map does not preserve the symplectic symmetry and does not have a constant 
of motion. To recover the symplectic symmetry we regularize the naive RG map by noting 
that the coefficient in the square brackets, multiplying A k (n) can be exponentiated: 



where 



u k 



AJn, 



A k (n + 1) = A k (n)exp iuJ A k (n) ) 



\ pi 1 ) \2 / a 

— h e — -Q ; cot cu fc + 2 smu k ^ Q ; Z> 



'(AT+l)^(AT) 





4iv\ 


A fc (n) 





JV=1 



(3.28) 



. (3.29) 



It is clear now that the regularized RG map ( |3.28| ) possesses the obvious integral of motion: 

Jk 



A k (n + 1) = A k (n) 



(3.30) 



It is worthwhile to note that the secular coefficients of the (2iV + l)-st harmonic e *( 2Ar +i)w fc n 
can be summed up to give a renormalized coefficient, which expressed in terms of A k (n) does 
not contain secular terms. 
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Proceeding in the same way as above, we can write the canonical conjugate momentum 
p n in the form 

p n = %B k {n)e lu)kU + c.c. + higher harmonics, (3.31) 

where 

B k (n + 1) = B k (n)exp iu k (^A k (n)^ . (3.32) 

Using now the relation (|3.1j) between the canonical conjugate variables (q,p) we can express 
the renormalized amplitude B k (n) in terms of A k (n) as 

Bk (n) = i- : C -^A k (n). (3.33) 

sin u k 



Neglecting higher harmonics and iterating Eqs. ( |3.28| ) and ( |3.32| ) we can write the renor- 
malized solution of the beam-beam map (|2.17|) 



q n = y2j£cos^ fc (J fc ;n), (3.34) 

p n = a k (J k )^2T k cos i(j k (J k ;n) - (3 k ( J fc )y / 2^4sin if) k ( J k ; n), (3.35) 

where 

ipk( Jk) n) = [uj k + ui k { J k )\n + 4> k , (3.36) 

cosu k - cos [u k + uj k ( J k )} afi\ sin[u k + u k (J k )} 

a k {J k ) = : , Pk{Jk) = : • (3.37) 

sin u k sin u k 

It is easy to see that the integral of motion J k has the form of a generalized Courant-Snyder 
invariant and can be written as 

2J t = f + (3.38) 

It is important to emphasize that Eq. ( p.38[ ) comprises a transcendental equation for the 
invariant J k as a function of the canonical variables (q,p), since the coefficients a k and j3 k 
depend on J k themselves. 

4 The Invariant Phase Space Density 

If an integral of motion J k of the beam-beam map ( p. 171 ) exists, it can be proved that the 
invariant phase space density f k (q,p) [which is a solution to the Vlasov equation ( |2.4|) ] is 
a generic function of J k , that is 

f ( k \q,p) = F k (J k ) (k = 1,2). (4.1) 
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Here Fk(z) is a generic function of its argument. Since the integral of motion J k is a functional 
of the invariant density of the opposing beam f^_ k (q,p), Eq. (|4.1|) comprises a coupled system 
of nonlinear integral equations for the invariant densities of the two counter-propagating 
beams. Let us find the integral of motion [see Eq. (|3.38|) 1 up to first order in the perturbation 
parameter e. We have 



Jfc — J() — 



\ k C% (J ) 



p cot Uk + pq), 



(4.2) 



where 



(4.3) 



The Fourier image of the beam-beam potential 



oo oo 



V k (X) = ~ J dq'J dp' f^ k (q',p') cos \q', 



— oo — oo 



(4.4) 



obtained by solving the Poisson equation (|2.2|) is next substituted into the corresponding 
expression [see Eq. ( |3.15|) 1 for the coefficient C^\Jq). Taking into account the recursion 
relation ( 3.19|) as well as the identities [10 



dxj'o(x) cos ax 



[0 < a < 1], 



(4.5) 



dxj 2 (x) cos ax 



2a 2 - 1 
VT^a 2 



[0 < a < 1] 



(4.6) 



we obtain 



oo \J2 Jq 



CPW = jJ d P'J dJf&tfrfWUo-q*- 



-oo 



(4.7) 



Thus, we finally arrive at the system of integral equations for the invariant phase space 
densities fj. (q,p) 



fi I \q,p) = C l F l 



oo \/2Jo 

J - (p 2 cot u)i + pqj J dp J dq'fP (q , p) ^2J - q' 2 



oo 



(4i 



A I \q,p) = C 2 F 2 



4 \ 00 ^ 2J ° 

Jo - -j 1 - (f cot u 2 + pq) J dp J dq'fP (q, p) ^2J - q' 2 



-oo 



(4.9) 
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where 



Oh 



oo oo 



-I -1 



dp / dqF k (q,p) 



(4.10) 



It is instructive to calculate the first order nonlinear incoherent beam-beam tune shift. 
According to Eq. ( |3.29| ) we have 

4\j k ) = ^Jdp' J djf£ k tf,i/)y/2J k -q*. (4.11) 

k -oo 

Since we are interested in the first order 0(X k ) contribution, we substitute in Eq. (|4.11| ) the 
unperturbed phase space density 

/f>(?,p) = ^p(-4±A (4.12) 



2<r 



Simple manipulations yield 

u k \ J k) 



: exp 



Jh 



dip(l + cos -0) exp 



Jk 



(4.13) 



Taking into account the integral representation of the modified Bessel function T n (z) (see 
e.g. Ref. 0]) 



J n (^) = _ / drcos(nr)e zcosr , 

7T J 




we obtain 



X k V2n 

0"3-A; 



■ exp 



Jk 



In 



Jk 



M-k. 



Jk 



M-kji 



(4.14) 



(4.15) 



A similar expression for the incoherent beam-beam tune shift was obtained in 



5 Concluding Remarks 

As a result of the investigation performed we have obtained a regularized symplectic RG 
beam-beam map, that describes correctly the long-time asymptotic behavior of the original 
system. It has been shown that the regularized RG map possesses an integral of motion, 
which can be computed to any desired order. The invariant phase space density (stationary 
distribution function) has been constructed as a generic function of the integral of motion, 
and a coupled set of nonlinear integral equations for the distributions of the two colliding 
beams has been derived. Based on the explicit form of the regularized RG map, the incoher- 
ent beam-beam tune shift has been computed to first order in the beam-beam parameter. 

It is worthwhile to note that the method presented here is also applicable to study the 
four- dimensional symplectic beam-beam map, governing the dynamics of counter-propagating 
beams in the plane transverse to the particle orbit. 
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